Analysis of the inelastic contribution to atomic spin excitation 

spectroscopy 

Bruno Chilian, Alexander A. Khajetoorians, Jens Wiebej*] and Roland Wiesendanger 

Institute of Applied Physics, Hamburg University, 
Jungiusstrasse 11, D-20355 Hamburg, Germany 

Abstract 

We investigate the dependence of the inelastic spin excitation contribution to scanning tunneling 
spectroscopy recorded above Fe adsorbates on the InSb(llO) surface at subkelvin temperatures, on 
bulk doping and tunnel junction resistance. To explain our observations, we show how the inelastic 
contribution depends on the parameters describing the excitation mechanism and tunnel conditions, 
in the framework of a recently developed model. We conclude that in this particular system of 
an adsorbate which is strongly relaxed into the substrate, the tip-sample distance dependent by- 
tunneling has to be taken into account in order to explain the observed variations in the inelastic 
tunnel spectra. 

PACS numbers: 68.37.Ef,72.25.-b,72.10.-d,71.70.-d 
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In the quest for nanometer scale spintronic devices, the scanning tunnehng microscope 
(STM) is a valuable tool due to its unrivaled ability to locally probe and manipulate the mag- 
netic structure of the sample system with atomic resolution. Most notably, inelastic scan- 
ning tunneling spectroscopy^"^ (ISTS) and spin-polarized scanning tunneling spectroscopy^ 
(SP-STS) experiments have contributed decisively to the understanding of single magnetic 
adsorbates on surfaces and their interactions^!^. 

Understanding the mechanism governing ISTS experiments has been the purpose of nu- 
merous theoretical advances^^-ii, which relate the experimentally accessible excitation spectra 
to the unobservable quantum processes that occur during tunneling. The descriptions range 
from a generalized Anderson model^i^ in a co-tunneling picture^^ commonly used in trans- 
port theoryi^ to approaches that treat the excitation event in a sudden approximation^!^, 
with a formalism similar to scattering theory. All of these models accurately reproduce 
the measured intensities of a given measurement, opening up the possibility of probing the 
hidden spin-excitation mechanism by comparison of model and experiment. 

In this work, we analyze ISTS and magnetic field dependent spin-split Landau level 
spectroscopy of single Fe adsorbates on the InSb(llO) surface containing a two-dimensional 
electron system (2DES)^'^ for differently doped substrates and for different junction resis- 
tances. We observe a significant variation in the inelastic contribution to our spectra with 
the doping concentration and junction resistance. To explain our observations, we utilize the 
model which treats the excitation events in a sudden approximation, where the formation 
and decay of an intermediate total spin state of tunneling electron and adsorbate governs 
the excitation probabilities^. We show for an arbitrary spin, how the inelastic signal per- 
centage Jinei and the 2DES Landau level asymmetry All depend on the model parameters 
that describe the excitation mechanism and tunnel conditions. By combining the measured 
Jinei and All, it is then in principle possible to extract the model parameters. We discuss 
the various physical effects that determine /inei and All and conclude that in our case, 
distance-dependent by-tunneling is most likely the origin of the variations in Jinei- 

All measurements were performed in an ultra-high vacuum (UHV) STM facility with a 
magnetic field of up to 12 T perpendicular to the sample surface (out-of plane) operated at 
a temperature of 300 mK—. We used electrochemically etched W tips, which were flashed 
in UHV to a temperature of approx. 2000 K. The n-doped InSb crystals with doping 
concentrations in the range of 0.5 — 2 ■ 10^^ cm"^ were cleaved in UHV to expose the clean, 
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atomically flat (110) surface^, onto which single Fe atoms were deposited at low temperature 
(T < 25 K)^. To obtain the dJ/d\^-spectra, the feedback circuit was switched off at a 
stabilization current /stab and sample voltage Vstah of 4 — 20 mV for the ISTS and 100 mV 
for the 2DES/Landau level spectra. The sample voltage V was ramped while recording the 
differential conductance signal via a lock-in technique with a modulation voltage Knod (rms 
value) added to V {f = 828 Hz). 

A 2DES is induced at the surface by the adsorbate-induced downwards band bending^. 
We varied the band bending by choosing samples with two different doping levels, one 
higher doped with a carrier concentration (at 77K) of 2 ■ 10^^ cm~^ and one lower doped 
with 6.5 ■ 10^^ cm"'^. The spectroscopic signature of the 2DES is a step in the differential 
conductance with an onset in the bulk band gap [— 240mV to OmV] which indicates the 
energy of the lowest subband Ei of the 2DES, as is shown in Fig. [l]^a,b). From the difference 
in the lowest subband energy, which is Ei ^ —30 meV for the high doped, and Ei ^ —55 meV 
for the low doped sample, the difference in band bending Abb ~ 25 meV between the two 
samples was estimated. Fig. [1] (c) / (d) shows a schematic sketch of the tunneling conditions 
for weak/strong downwards band bending. Schematic sketches of the minority (pmin) and 
majority (pmaj) local vacuum densities of states (LDOS) above the Fe adsorbate, i.e. at 
the location of the tip, are included. Figures l(a,b) also show d//dl^ curves measured in a 
magnetic fleld of -B ~ 5 T revealing spin-split Landau levels (LL) as discussed later. 

The results obtained from the ISTS measurements for the high and low doped samples 
are shown in Fig. |2] (a,b). For both samples, the spectra reveal two steps above and below 
the Fermi energy Ep {V = OV), which are symmetric in energy with respect to Ep. These 
steps are due to the inelastic excitations of the adsorbate spin by the tunneling electrons^. 
However, the intensity of the steps strongly differs between the two samples, i.e. the inelastic 
contribution /jnei to the spectra differs. The spectra were normalized by dividing them by 
a spectrum taken on a nearby substrate location. Jinei was extracted from the normalized 
spectra by fitting a thermally broadened symmetric double step function and then taking 
the difference of the function above {V > 5 mV) and below (V^ = mV) all excitations [gray 
horizontal line in Fig. |2] (a,b)]. For the higher doped sample, measured with a junction 
resistance R = = 30 Mf2, we obtained /inei = 31 % ± 4 % (standard deviation) and for 

-'stab 

the lower doped sample, measured with R = 50 MQ, Jinei = 16 % ± 4 %, when analyzing a 
total of 24 adsorbates. To be able to separate the influence of doping and junction resistance, 
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we recorded /mei as a function of R on the higher doped sample (Fig. [2]^c)). We observe that 
at a junction resistance = 50 MTL on the higher doped sample, /inei ~ 20 — 25 %, which is 
only slightly larger than the Jinei measured on the lower doped sample at the same junction 
resistance. Considering the error margins, this implies that changing the doping may have 
no independent effect and the change in Jinei is entirely due to the higher R used on the 
lower doped sample. 

In ref. 17|], a technique complementary to ISTS for probing the magnetic properties of an 
adsorbate on the InSb(llO) surface was introduced. The asymmetry of the spin-split Landau 
level (LL) peaks of the Fe-induced 2DES was shown to reflect the adsorbate's magnetization 
as a function of magnetic field applied perpendicular to the surface. Fig. [2] (a),(b) shows 
typical d//dV^ spectra, utilizing a W tip with negligible spin polarization, taken above the 
Fe adsorbate in a magnetic field for the higher and lower doped samples. Both spectra 
show pairs of peaks with intensities d//dV^ and d//dV^ in the entire voltage range above Ei 
which are due to the spin-split LLs as marked by the corresponding arrows indicating the 
direction of the electron spin in each LL. Across the whole energy range, the spin-split peaks 
of each LL exhibit a large asymmetry All = dj/dy^+d//dvj[ (approx. 30 % for the lowest LL at 
i? = 5 T), with the lower energy spin-split peak having the lower intensity for both samples. 
As shown in ref. [l7|, this means that Fe majority spins are preferentially transmitted in the 
entire energy range from — 50meV to +100 meV. 

To understand the variations in the measured inelastic contribution to the ISTS spectra, 
we adopted the generalized version of the model used in ref. jol, which we describe in the 
following. The spin-dependent part of the Hamiltonian H for a magnetic adsorbate on a 
surface can be approximated hj H = —giiB^-S+DS^+E^S'^ — Sy) where S = (Sx, Sy, Sz) 
is the vector containing the operators of the adsorbate spin, B is the magnetic field and D 
and E are the magnetic anisotropy energies resulting from the spin-orbit interaction of 
the adsorbate spin with its surroundings. In the experimental situation discussed here, the 
magnetic field always points along the y-axis (out of plane), i.e. B-S = BySy. Consequently, 
we chose this the spin quantization axis. 

In ISTS, steps in the differential conductance appear due to the opening of additional 
tunneling channels when the tunneling electron's energy equals the energy difference between 
the magnetic ground state and an excited spin-state of the adsorbate. Consequently, D and 
E can be determined directly from the measured step energies. The step heights, on the other 



hand, depend on the excitation mechanism. Extracting the magnetic information encoded 
in the step heights thus requires a quantum mechanical model that links the underlying 
excitation mechanism to the experimentally observable quantities^. 

In the model proposed in refs. js] andlsl, the tunneling electron forms an intermediate 
total spin state with the magnetic adsorbate's spin S during the excitation process. This 
state has spin J = S* ± |. The adsorbate states are elements of the spin-S* state space 
Tis and the tunneling electron spin states are elements of the spin-| state space 'Hi/2- Let 
{(f) j \j = 1, . . . , 25*+ 1} be an orthonormal basis of Hs, where the (pj are eigenstates of H. In 
the usual notation, we define (0, m| = (0| ® (m| G Tis ®Ki/2 for (0| G Tis and (m| G 7/i/2- 

The term 

X.,m,f,n.' = I (0„ m| c a+ + (1 - c) e*'^ a_ m') |2, (1) 

with c G [0,1] and (/? G [— 7r,+7r], describes the relative transition intensity of a tunneling 
event in which a tunneling electron in the initial state \m) interacts with the magnetic 
adsorbate in the initial state After the tunneling event, the adsorbate is in the final 
state 10/) and the electron is in the final state |m'). These transitions are governed by the 
operators 0± = XIa/j I'S' ± V2, Mj){S ± 1/2, Mj\, which describe the tunneling event as the 
formation of an intermediate state | J, Mj) with total spin J = ^±1/2 and magnetic quantum 
number Mj. The strength of the transition is proportional to the product of overlaps of the 
initial states with the intermediate state and the intermediate states with the final state. The 
summation over the different unobservable intermediate states is performed before taking 
the absolute square, so the different tunnel paths interfere. Note, that this model allows for a 
situation, where the electron transmission is governed by a mixture of majority (J = 5*+ 1/2) 
and minority {J = S — 1/2) spin orbitals in the LDOS of the Fe adsorbate. The parameter 
c governs the relative strength and (f the relative phase of the contribution of these two spin 
channels. 

Before the model is analyzed in detail, we compare to other approaches reported in the 
literature. In refs. and |2l|], a similar model was used, where the relative transition 
intensities are given by the expression |(0j,m|S ■ s + u ■ l|0/,m')p where s = {(Jx,cry,az) 
is the vector containing the electron spin Pauli matrices, 1 is the identity operator. From 
the properties of 0± given below, it is possible to derive the relation A+O^ + A^O^ oc 
S • s + ^^('^^^^^)~ ^ ; which shows that the transition operators defining the two models can 
be mapped onto each other, with the transformation from {A^, A_) to u given by the second 



term on the right hand side. 

We proceed by analyzing /jnei and All as predicted by these models. At low temperature, 
only the ground state will be occupied if we neglect artificial pumping of the spin states 



by the tunnel current as justified in ref. 17|. If the ground state is non-degenerate, the 



elastic contribution to the transmission is proportional to Xgi = '^m-^9,m,g,m- This term 
describes all tunnel events in which the adsorbate's spin remains in the ground state during 
the collision, i.e. there is no excitation. Correspondingly, the transmission at energies higher 
than all excitations is proportional to Z = J^m'^fm' -^g,mj,m', which describes all allowed 
events where the adsorbate starts in the ground state. The central quantity of interest in 
this work is the fraction of the differential conductance that is due to inelastic tunneling 
events, Jinei = ^ ~ This expression can be generalized to the case of a degenerate 
ground state by summing over the different ground states in Xei and Z. For simplicity, 
we will assume a non-degenerate ground state in the following, as this is the case in our 
experimental situatioi>i^. 

Using the fact that the intermediate states \J,Mj) form a complete and orthonormal 
basis set of Ks ® K1/2, one can show that 6+ + = 1, O+O. = 0, = d±, and 
S ■ s = \SO+ — ^{S + 1)0-. The expectation value of the spin scalar product (S ■ s)j^m = 
{(pi, m|S-s|0j, m) is an important quantity in this model as will become clear below. Because 
all terms in H except the Zeeman energy are quadratic in the adsorbate spin operators, it 
is zero at B = 0. It reaches its maximum magnitude max|(S ■ s)j^m| = ^/'^S when the 
eigenstates of H are simultaneous eigenstates of Sy, and it is evaluated for an eigenstate 
with the magnetic quantum number M = ±S. It satisfies (S ■ s)^ = — (S ■ s)j _i . 

Using the above relations, one can show that Jind has the following properties for all 
D, E, B: (i) phase symmetry /inei(c, (p) = /inei(c, —^p)- Therefore we limit all further investi- 
gations to the domain (f G [0, vr]. (ii) A global minimum Jind = occurs at c = 0.5, (f = 0, 

(iii) a global maximum /inei = s{s+i) ^' occurs at c = ^/2S+i,(p = n. For B = 0, this 

is maximal, Jinei = 1, and if the adsorbate's spin is fully aligned to the magnetic field, this 
maximum reaches its smallest possible value ^s+i. (iv) There are no further extrema and 
/inei increases with ip, except at the boundaries c = and c = 1, where /inei is independent 

of ip, Jinei (c = 0) = — 2Sis+^)''^ ' ^''^^''^^ = 1) = 2is+i)(s+ff ' ^hcsc results imply that 
for B = and at the maximum possible phase difference ip = n, if the intermediate spin 
channels contribute with just the right mixing c = ^/2S+i between majority and minority 



spin orbitals, the elastic transmission is completely suppressed and the excitation efficiency 
is 100 %. On the other hand, for all magnetic fields, if the channels contribute equally with 
no phase difference, all excitation channels cancel out and only elastic tunnel events occur. 



The behavior of /jnei as a function of c and ip for the experimentally relevant S = 1 and 
the anisotropy energies that were found for Fe on InSb(llO)^^ at B=0 is shown in Fig. [3t^a). 
The global maximum is located at c = 1/3, = vr, and /jnei takes on all values in [0, 1]. Note 
that every contour of constant Jinei intersects at least one of the boundaries = 0, = vr, 
so any value of Jinci can be produced by real coefficients of 0±. 

It is possible to show that the ratios of the inelastic steps ^, Xg^m,h,m' / Xg,m,f2,m' 
for fi,f27^g are independent of c and ip, so these do not provide an independent source of 
information about the excitation mechanism. We therefore conclude that knowledge of Jinei 
at one magnetic field alone only implies which contour on the c-y^-plane describes the tunnel 
conditions. To further narrow down the range of parameters, a complementary source of 
information is required. What is lacking is an observable that refiects the evolution of the 
adsorbate state in a magnetic field. 

The spin-split LLs shown in Fig. [II^a,b) provide a means of selecting the initial spin state of 
the tunneling electron, by setting the (negative) sample bias voltage to either the position of 
the spin-up or the spin-down peak. If the completely spin polarized tunnel current from one 
of these peaks is passed through a magnetic adsorbate on its way to the tip, it may interact 
with the local moment, creating spin excitations of the adsorbate and fiipping the tunneling 
electron spin as a result. If the local moment is (partly) aligned by the magnetic field, this 
creates asymmetric tunnel conditions for the different tunnel electron spin orientations. As a 



result, the observed intensity of the two spin split LL peaks will be different. In ref. [17|], the 
model ([Tj) was adapted to describe the LL asymmetry ^LLfor the case S = 1. Using a similar 
analysis as for Jinei, it is possible to derive a rather general expression for All- Let p+ be 
the relative peak height that corresponds to tunnel electrons that have spin up in the initial 
state and let p- correspond to spin down. These quantities are given by the total transition 
intensity for the adsorbate in the ground state and the electron in the up (down) state to be 
transmitted into any other configuration: p± = X]/m' I (0c/5 =1=^1 CO4. + 6**^(1 — c)^-!^/, m')p, 
since the energy of the lowest LL is much greater than all excitation energies, i.e. all channels 



are open. With this, we get the LL asymmetry 

-2(2c-l) 

"^"^ P-+P+ c^{2S + l) + S{l-2cy ^^•+^' 

This is proportional to the projection of the adsorbate's spin on the magnetic field axis 
{4>g\Sy\4>g) , i.e. to the magnetization of the adsorbate. All undergoes a sign change at 
c = 0.5, so the sign of All determines whether contours in the upper or lower half of 
Fig. [3](a) have to be considered. 

Fig. |3](b) shows the predicted asymmetry curves for the experimentally relevant param- 
eters of Fe adsorbates on InSb(llO)^^ and several values of c. We note that in principle, if 
All were known experimentally with high precision, c could be determined precisely from 
Fig. El^b) and then, using the the measured Jinei and Fig. [31(a), ip would be known, too. This 
demonstrates the complementarity of Jinei and All- However, for c close to either 1 or 0, the 
asymmetry curves are hard to distinguish experimentally, so it is unpractical to determine 
c with high precision. From the positivity of our measured All, we deduce that c > 0.5, i.e. 
the S + ^ total spin channel dominates the tunneling. In other words, the majority spin or- 
bitals dominate the Fe vacuum LDOS in the whole energy window. From the measured Jinei 
for the higher/lower doped sample at i? = 30Mn/50Mfi, we conclude that the parameters 
describing our system lie in the area highlighted in green/yellow in Fig. [3]^a). This result 
shows how in principle the measured Jinei and All can be used to extract the parameters c 
and (f. 

In the following, we discuss the results. The value of the parameters c and (p are de- 
termined by the details of the excitation process, namely whether tunneling electron and 
adsorbate spin preferentially couple parallel or antiparallel to form the intermediate total 
spin state and, if both channels are present, with which relative phase they contribute. Par- 
allel coupling of the adsorbate and electron spins is facilitated if the vacuum LDOS above 
the adsorbate is dominated by majority spin orbitals and antiparallel coupling is promoted 
by a minority dominated vacuum LDOS. In this way, c refiects a crucial magnetic property 
of the adsorbate's electronic structure. Additionally, the formalism can be used to describe 
so called by-tunneling, where the tunnel electron does not couple to the adsorbate at all 
but bypasses the magnetic adsorbate orbitals, leaving the adsorbate and electron in their 
initial spin states. This can be achieved by adding a multiple of the identity operator 1: 
cO^ + (1 — c) e**^ 0_ + b - 1. But since i = (9+ + (!)_, this additional term can be absorbed 
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into the coefficients of the 0± and the shape of the transition operator remains unchanged, 
a(c'(9++e*'^'(l — c')(!)_), where the scale factor a has no effect on the relative transition inten- 
sities. This implies that the effect of by-tunneling on /jnei and All cannot be distinguished 
from spin channel mixing. However, as we will discuss next, in our case, by-tunneling is 
better suited to explain the experimental observations. 

In principle, the variation of band bending or junction resistance can influence /inei in the 
following ways: (i) As illustrated in Fig. [T]^c) / (d), the spin polarization of the vacuum LDOS 
at the Fermi energy may be energy dependent, resulting in changes in the mixing parameter 
c as Ep varies relative to the adsorbate LDOS. Though an intriguing possibility for tuning 
the excitation efficiency, this is unlikely to be the case in our experiment, since the difference 
in band bending (Abb = 25meV) is too small to produce a. n/A change in (y9 or a change 
of approx. 20 % in the spin-polarization, which are the values that could account for the 
experimentally observed change in Jinei (see Fig.3(a)). (ii) Since the 2DES density is larger 
for the lower doping and tunneling into the substrate is presumably dominated by the 2DES, 
the tip-sample separation will be larger for lower doping. Likewise, increasing the junction 
resistance R will increase the tip-sample separation z. If the vacuum spin polarization varies 
with distance from the surface, this may account for changing c. However, it seems unlikely 
that the vacuum spin polarization changes sufficiently on the distance scale of a few tens 
of pm^^, which is a rough estimate of the relevant change in z by i? oc exp(2z/A). A more 
plausible explanation is that (iii) the by-tunneling contribution may be strongly distance 
dependent, because the Fe adsorption site lies below the topmost In and Sb layer—. Our 
experiment thus shows, that by-tunneling is getting stronger for larger tip-sample separation, 
which seems plausible. 

Our method of the determination of c and (p that describe the excitation process will be 
most effective in a setting where by-tunneling is known to be minimal, so that the pure spin 
channel mixing can be observed. This is expected to be the case for adsorption geometries 
where the magnetic adsorbate is not strongly relaxed into the surface. We note that the 
measured Jinei = 31 % ± 4 % is compatible with the value /inei(c = 1) = = | for S* = 1, 
B = 0. This on the one hand justifies the choice of c = 1 in ref. 17| and on the other hand 
shows that for low junction resistances on the higher doped sample, there is practically no 
by-tunneling. 

In conclusion, we have shown how the inelastic contribution to ISTS spectra li^ei and the 



Landau level asymmetry All can be used to investigate the mechanism of the spin excita- 
tion process. Our results motivate further investigations into the possibility of tuning the 
excitation efficiency. In particular, it would be desirable to obtain an independent measure 
of the by-tunneling with spatial and energetic resolution, as this cannot be distinguished 
from basic elastic tunneling through the magnetic adsorbate at this stage. 
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FIG. 1. (a) Typical 2DES and LL spectra for the high doped sample, exhibiting weak band 
bending. The 2DES (black) and the LL (red) spectra, were taken in zero magnetic field, and in 
a magnetic field By perpendicular to the surface, respectively, (/stab = 200 pA, Kitab = 100 mV, 
Vmod = ImV). (b) The same for the low doped sample exhibiting strong band bending, {1^^^^ = 
200 pA, /g^fb = 150 pA, Fstab = 100 mV, V^od = 0.5 mV). The 2DES spectrum was multiplied by 
3/4 to compensate for the different /stab- The spin states of the spin-split LLs are indicated by 
arrows. (c)/(d) Schematic representation of the band bending for high/low doping. The conduction 
band (CB) and valence band (VB) edges are energetically lower at the surface (vertical line) than in 
the bulk, inducing a surface 2DES. An artificial distribution of the majority and minority vacuum 
LDOS of the Fe adsorbate is sketched. 
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FIG. 2. (a)/(b) Typical ISTS spectra of isolated Fe adsorbates taken at By = (T = 0.3 K). The 
spectra were normalized by division by a spectrum taken at a nearby substrate location. Broad red 
curves represent the best fit to a thermally broadened symmetric double-step function and gray 
horizontal lines mark the inelastic contribution /jnei to the total differential conductance (see text), 
(a) Higher doped sample (/stab = 600 pA, Kitab = 20 mV, Vmod = 40 /uV). (b) Lower doped sample 
(/stab = 200 pA, Fstab = 10 mV, Vmod = 40 ^V). (c) /inei as a function of junction resistance R on 
the high doped sample. The error bars were computed from a 90 %-confidence interval to the fit 
parameters. 




FIG. 3. Model predictions for an adsorbate with spin 5=1 and anisotropy energies D = —1.4 meV, 
E = 0.22 meV. (a) Contour plot of /jnei as a function of the model parameters c and ip. Shaded 
areas highlight the experimentally relevant range of parameters for high/low doping in green/yellow 
(b) Landau level asymmetry ^ll as a function of magnetic field By along the y-axis for various 
values of the mixing parameter c. 
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